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Abstract 

Two interacting electrons in a harmonic oscillator potential under 
the influence of a perpendicular homogeneous magnetic field are con- 
sidered. Analytic expressions are obtained for the energy spectrum 
of the two- and three-dimensional cases. Exact conditions for phase 
transitions due to the electron-electron interaction in a quantum dot 
as a function of the dot size and magnetic field are calculated. 

PACS numbers: 73.20 Dx, 73.23.Ps 



1 Introduction 

In recent years considerable experimental and theoretical interest has been 
focused on the study of artificially structured atoms in semiconductors, such 
as quantum dots, where few electrons are confined in all three directions (see 
for review |I|, ^, ^). In experimentally realized quantum dots, the extension 
in the x — y— plane is much larger than in the z direction. Assuming that 
the z-extension could be effectively considered zero, the electronic proper- 
ties in these nanostructures have successfully been described (see |jl|, Q and 
references therein) within the model of the single-electron motion in the two- 
dimensional harmonic oscillator potential in the presence of a magnetic field 
[^]. Based on a numerical solution of the Coulomb interaction between elec- 
trons, a complex ground state behavior (singlet— ^triplet state transitions) as 
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a function of a magnetic field has been predicted [|, ^ (see also discussion in 
[0]). Remarkably, these ground state transitions for N=2 have been observed 
experimentally 

In present paper we consider an analytically solvable model of two elec- 
trons in a quantum dot. The confinement potential is approximated by 
a harmonic oscillator potential and the problem of the Coulomb interac- 
tion is treated exactly. Though this case represents the simplest nontrivial 
problem with regard to the electron number, mainly the ground states of a 
two-dimensional quantum dot were analysed either quantitatively or an- 
alytically with some approximations [§. The role of third dimension {z) is 
also investigated and results on the analysis of the ground as well as excited 
states are presented. 

2 Model: General Remarks 

Our analysis is based on the oscillator representation method (ORM) devel- 
oped in IQ]. The ORM arisen from the ideas and methods of quantum field 
theory has been proposed to calculate the binding energy of different systems 
with fairly arbitrary potentials described by the Schrodinger equation [|10|] . 
Here, for completeness we present briefly the main ideas of the ORM. 

For any potential admitting the existence of a system bound state there 
is always a transformation of the variables that leads to a Gaussian asymp- 
totic form for the wave function at large distances. However, the asymptotic 
behavior of the Coulomb wave functions for large distances does not coin- 
cide with this behavior. Therefore, we have to modify the variables in the 
original Schrodinger equation so that the modified equation should have solu- 
tions with the Gaussian asymptotic behavior. In the Coulomb systems, this 
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modification is performed by going over to the four-dimensional space, where 
the wave function of the Coulomb system becomes the oscillator one. In an 



early paper |TT|, Schrodinger has noted the existence of such a transformation 
which transforms the three-dimensional Coulomb system into the oscillator 
one in the four- dimensional space. The explicit form of this transformation 



has been found in [|T^ and used to solve the classical Kepler problem. 

In the next stage, it is necessary to represent the canonical variables 
(coordinate and momentum) of the Hamiltonian through the creation and 
annihilation operators and a. From the Hamiltonian the pure oscillator 
part with some, yet unknown, frequency uj is extracted, i.e. H => Hq + Hj = 
ujo^a + higher order terms. The remaining part, i.e. the interaction 
Hamiltonian Hj, is represented in terms of normal products over a"*" and a. 
In addition, it is required that the interaction Hamiltonian does not contain 
terms quadratic in the canonical variables. This condition is equivalent to 
the equation 

^ = (1) 

which determines cu, the oscillator frequency, in the ORM and is called the 
oscillator representation condition (ORC) §]. Similar ideas are used in the 
Hartee-Fock-Bogoliubov theory to describe different correlations between nu- 
cleons moving in average nuclear potential (see for review |p!3|). 



Since we change our space variable r and a magnetic quantum number m 
will be absorbed by the dimension parameter of the auxiliary space (0 and 
see below), the calculation of the wave function 'il){r) would be equivalent 
to the calculation of the ground state function of a modified Hamiltonian in 
another dimension. As it has been mentioned above, the wave functions in 
this auxiliary space should have the oscillator Gaussian asymptotic behavior 
at large distances. This property is quite natural for our purposes due to 
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physical conditions for the confined electron gas in a quantum dot. The 
effective confining potential (oscillator) prevents the tendency caused by the 
Coulomb forces to allow electrons to escape and, therefore, it should dominate 
in the phenomena. 

The model is described by the Hamiltonian 



H 



E 



2m' 



e -> .0 
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+ 



spin 



Arceeo |fi - 

where Hspin = g{si + S2)B. Here m* is the effective electron mass. Below, 
we use the units (e = c = 1). For the perpendicular magnetic field 
we choose the gauge described by the vector A = [B x r\/2 = ^B{—y,x,0). 
Introducing the relative and center-of-mass coordinates 



fs-ri R=-{fi + f2) 



(3) 



the Hamiltonian, Eq.(|^), can be separated into the center-of-mass (CM) and 
relative- motion (RM) terms as (see also P, ||) 



H 



Ir^ -<-|2 ^ ^ 2 2x k^/hujQ 

h 



Q 



Pq + Aq + [ujq'p'q + u'q^QI 



(4) 
(5) 

(6) 



where ujq = 2uJo, uq^ = 2uj^, lo^ = ^Uq, uj^^ = ^uj^, Aq = A{qi) + A{q2), 
Aq = ^ (^A{q2) — A{qi)^ and A{q) = :^[B x q\. Here we have introduced the 
variables q 



-r,Q 



lengths: the effective radius a 



R, p = y/x^^+y^ and defined the characteristic 



and the oscillator 



length Iq = (h/m*uJoy^'^. These units allow one to define the dimensionless 
dot size k = lo/a* 0. 

The separability and the conservation of the angular momentum lead to 
a natural ansatz for the eigenfunction of the Hamiltonian, Eq.(|^), 

^ = HmQ)xisuS2), (7) 



where the wave functions i/j{q) and (f){Q) are [[14 



'ip{a) = —=i)^{pa, z) (8) 



/27r 

and the eigenvalues have the form 

E = 2er + -E]y^M + Espin- (9) 

Here and E^^m are the eigenvalues of the Hamiltonians Hg and Hq, re- 
spectively. According to the Pauli principle, if the spatial part of the total 
wave function is symmetric (antisymmetric) with respect to the inversion 
r — r, X must be the singlet (triplet) spin state. 

We now concentrate our analysis on the relative motion Hamiltonian Hj.. 

3 Coulomb problem 

Due to the axial symmetry of the problem, the Schrodinger equation with 
the Hamiltonian, Eq.(^, can be written in the form 



d'^ Id w? 

dpi Pq dpq PI 



1 r/2 f)'^ 

I ^ /O 2 2 I 2 2\ 
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+ :rr^=^^ra{pq, Qz) = f^m^m(Pq, Qz) (lO) 



2a/p + q. 
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Here 



(11) 



where t = Uc/ujq, uic 



is the cyclotron frequency, and 




(12) 



where m is a magnetic quantum number. 

According to the ansatz of the ORM for the wave function, we have to 
change radial variables so as to obtain an oscillator asymptotic behavior for 
the wave functions of the transformed equation and then identify this equa- 
tion with the radial Schrodinger equation in a space with different dimension. 
In addition, our Hamiltonian contains a repulsive centrifugal term, and the 
wave function ipmiPq, Qz) niust decrease at small distances. Consequently, the 
transformation to the higher dimensional space is realized by 



The parameter d can be chosen to compensate completely the repulsion at 
small distances. The calculation of the wave function $m.(Pg, Qz) is equivalent 
to the calculation of the ground state wave functions in a space W^. The wave 
function, Eq.(|l3D, is a regular one at short distances. Therefore, our wave 
function possesses necessary properties at short and large distances, i.e., it 
is a Gaussian one as p — oo and goes to zero at p — » 0. 

According to the definition of the wave-function, Eq.(0), we can trans- 
form Eq. (p!OD into the equation 



(13) 




(14) 
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which can be identified with the equation in space with 



d = 2 + 2\m\ (15) 

One can see that the magnetic quantum number m does not enter into the 
Schrodinger equation, Eq.(0), in the exphcit form. It is absorbed by the 
"dimension" parameter d. This trick allows one effectively to avoid the prob- 
lem of calculation of excited states and to perform calculations of the ground 
state in the auxiliary space R. Therefore, Eq.(p!4D contains the oscillator 
with the coordinate Pq G and the other one with the coordinate & R^, 
respectively. 

Choosing different (fixed) values of lj^ allows to study the dependence of 
the results on the slab thickness. The condition Uz 3> ujq ensures that we 
have a genuine two-dimensional problem implying that no particles occupie 
a quantum mode in the 2;-direction (n^ = 0). From the analysis of the far- 
infrared frequencies in the three-dimensional oscillator potential it follows 
IH] that the frequency that just forbids occupation of a 2;-mode is given by 



u'^ > Lj±{y/4:N + 1 — 3)/2 with u± being the average of and Uy. 
3.1 Two-dimensional case 

Let us consider the case z = 0, i.e., the pure two-dimensional electron gas. 
For the case under consideration Eq.(|T^ can be written as 

H2d^m{pq) = f/m$m(Pg) (16) 

where 



1 

Horl = 

2d 2 



(i^ d — ld 
+ 



+ R + ^ (17) 



dpi Pq dpq\ 2 2p, 
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Here the wave function $m(pg) depends only on pq = Therefore, we 

can identify the operator 

d'^ d — 1 d _ ^ 
dpi pq dpq 

with the Laplacian Ap^ in auxihary space R'^ if this operator acts on a func- 
tion depending on the radius only. The wave function $m(pg) in Eq.(|r^) 
can be considered as a wave function of the ground state satisfying the 
Schrodinger equation 

H<^m{p,) = e{E)<^>M (19) 

where 

H = ^ + %m + ^ - (20) 

Taking into account Eq.(|l6]), the desired energy E is determined by the 
equation 

e{E) = (21) 

Let us express the canonical variables p and p through the creation and 
annihilation operators and a 



1 

(22) 




where is a new oscillator frequency which is defined below. The vacuum 
state |0 > is defined according to the standard rules 

<0|0>=1, aj|0>=0 (23) 
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Substituting the representation, Eqs.(0), into the definition of the Hamiho- 
nian, Eg. pOj) , after some transformations 



^ fi.5:a;a, +A(^ + ^^) (24) 



and 



Pq J J y/TT J y/TT 

— OO — CXD 

/^u; (^^ffll +J^I (^) ^"^(1+^^) : ^-^(^^"^ : j (25) 



-OO 



we obtain 

H = Ho + Hi + eo (26) 



where 



i/o = hujj^ajaj (27) 



^ dnl\ h ^ r(^^) 

Hi = ^y^khj (29) 



Here : ... : means a normal product, and we have introduced the notation ef = 
— 1 — X — According to the ORM f^, the interaction Hamihonian Hi 
does not contain terms quadratic in the canonical variables, i.e. proportional 



to : pi :. 
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The ORG requirement, Eq.(|I|), determines the oscillator frequency uj 
which is defined in the following way: 



UJ = X Vtr, 



(31) 



The quantity x is determined with the following equation which is derived 
from Eqs.dH), (0) 



X 



k r(i + |m|) 



-1 = 0. 



(32) 



^ (1 + 1)1/4 r(2 + |m|) 

It is clear that at zero Coulomb field {k = 0) x = Xq = 1 while for k 
Eq.(l32|) defines the effective dependence on the Coulomb interaction of the 
oscillator frequency {k « 1). Considering the quantity x expanded as 
a Taylor series in the variable k and keeping only first order terms x = 
Xq + kxi + ... , we obtain, according to Eqs.(|ll|), (pl|) and (|32|) 



t2 



t2 



UJ 



where 



UJa 



UJa 



[I + kXl)uJg\l 1 + J = ^^q\J 1 + — 



l/a* r(i + |m|)^ 



(33) 



(34) 



2v^(l + it2)i/4r(2 + |m|); 
When the Coloumb forces are absent {l/a* = 0) it follows that ujg = ujq and 

UJ = Vlq. 

According to the ORM, the quantum number n defines the radial excita- 
tion (see P), i.e., the highest oscillator states 



Cn 



> 

1/2 



(35) 



_22"n!r(f + n) 

Correspondingly, the energy spectrum with radial excitations is defined as 

e^''\U)=<n\H\n>=ai + a2 (36) 
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with 

ai = {- + 2n)huj+'^^-Um (37) 

n ( r(— )\ 

a2 = -V^k n\hi\n > 
Taking into account Eqs.(0), (|32D , ( ^61) from Eq. (|2l|) we obtain 

nm '^nm ' ^nm 



_ M) 

2 



^nm 



^t + (l + |m| + 2ny (1 + ^)1/2 



(38) 



£ 

nm 



2 2V2^ a' 



' r(i + |m|) , 

r(l + \m\) 



where the matrix < n\hj\n > is defined in Appendix A. In perturbation the- 
ory the effect of Coulomb forces is taken into account by the second term e^^. 
In our approach the main term depends on the Coulomb forces as well. 
In standard schemes this term corresponds to the noninteracting electrons 
moving in the external confining potential [Q. Here, within our model, the 
interaction modifies the external potential and results in the effective mean 
field potential of the relative motion. 

3.2 Three-dimensional case 

Despite that Uz ^ ujq, in real samples the effect of the third direction should 
be taken into account, and the prediction based on the pure two-dimensional 
case is expected to be modified. 

Taking into account the definition Eq.(|r^), Eq.(|T^) can be written in the 
following form: 

[{H2cl - Um) + hz + Kes] $m(Pg, Qz) = (39) 
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where 



1 Tf 




(41) 



'Pi + Qt 

Since the terms and give the main contribution to the total Hamiho- 
nian, the term hres is related to a dimension of the problem and can be con- 
sidered as a correction term. Let us introduce a transformation for the one- 
dimensional oscillator hz similar to the two-dimensional case (see Eg. (p2D ) 



After some transformation of the Hamiltonian, Eq.(pU|), we obtain 



Here 



H = Ho + eo + ^.J^okhi (43) 



Hq = hUq^A^A + huo^a^aj (44) 

[ dTil + r')-'/'{l+^r')-'/' (45) 



20F 

and hj consists of four terms (see Appendix B). Here u is defined by Eg. (pi]) 
and 7 = — << 1. Finally, applying the definition of the radial excitations 
(see Eg . (1361) ) for the three-dimensional case, Eg. (^31) , from the condition. 



Eg. (|2l|) , we obtain the following expression for the lowest energy level with 
nz = 

hxk 



rp rp2d , 



f)"'Q(7)+2ia±Ms„(,) 

4 f (1 + v^\j 



(46) 



where the guantities (^(7), SJypi) are defined in Appendix B. 
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4 Discussion 



The solution to the Hamihonian of the center-of-mass motion Hq is well 
known and the energy can be written as 



where N = 0,1,... and M = 0, ±1,... are radial and azimuthal quantum 
numbers, respectively. The spin of the two electrons leads to an additional 
Zeeman energy 



m is a magnetic quantum number corresponding to the relative motion and 
g* is an effective Lande factor. 

Summing Eqs.(|^), (|48|), (^) (or Eq.(^6D in the three-dimensional case, 
respectively) we are able to investigate different ground states as a function of 
the dot size k = lo/a* and relative strength of the magnetic field uJc/ujq. Since 
the center-of-mass quantum numbers A^, M and the quantum number m are 
conserved by the Coulomb interaction, the ground state has the quantum 
numbers N = 0, M = 0, n = 0. Comparing the energy with different m < 
we can define the ground state energy for a given dot size k at different 
strength of a magnetic field ujc/ ujq. 

In our calculations, we used the effective mass m* = 0.067me of typical 
quantum dots for GaAs. In Fig.l (a) the energy spectra without a contribu- 
tion of the Coulomb forces are presented. While without the Coulomb forces 
the ground state is always the state with m = 0, the Coulomb interaction 
(Fig. 1(b)) leads to a sequence of different ground states m = — 1, —2, ... which 
are an alternating sequence of singlet and triplet states. 



= 2huo (2iV + |M| + l)Wl + - + (2n, + -) 




+ -Mt 



(47) 



1 ?77* (i) 

Es = g'^isS, = T [1 - {-m 

4 rUe ujQ 



(48) 
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The main mechanism, which defines the optimum quantum number m 
of the ground state, is the interplay between the dot size and the strength 
of the magnetic field. This question has nicely been discussed for the two- 
dimensional quantum dots in |^ (see also f^). Similar behavior is observed 
for the radial excitations with n = 1, 2, ... 

If the third extension (z) is taken into account, the ground phase transi- 
tions are shifted to a higher magnetic field (see Fig. 2). Since the extension 
of the slab is inversely proportional to the confining frequency {uJz ~ the 
thicker the slab the larger value of the magnetic field is needed to observe 
the ground state transition m ml. This fact has to be taken into account 
in experiments. 

The singlet-triplet ground phase transition occurs when the following 
condition is fulfilled Eo^rn = -£'o,m-i("^ < 0). For a negative Lande factor 
the spin-splitting energy in a magnetic field will lower the energy of the 
spin Sz = +1 component of the triplet states. In particular, the relation 
i?o,m = -Eo,m-i = -So,m-2 (m odd) defines the point when the singlet phase 
ceases to exist Beyond this point we can observe phase transitions be- 
tween triplet states defined by the condition -Eo,m = -Eo,m-2 (m odd). There- 
fore, at strong magnetic field Uc » ujq, i.e. in the limit t ^ oo and x ^ 1, 
for singlet-triplet phase transitions m —>■ m — 1 we obtain 



lo 8r(2 + |m|) 



fstil) (49) 



a* 3r(i + |m|) 

/-(7) = 1 + of!w'"' , . + 0(f) (50) 



and for triplet-triplet phase transitions m — m — 2 (m odd) 

h _ §^^0^3/2 ^(3+ |m|) 4 ^ 
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7 + 6 


m 




2 (5 + 4 


m|)(3-2 


m 


)(l-2 


m ) 



fuh) = 1 - , Jo Tin . + (52) 



In these expressions the pure two-dimensional case is reahzed in the hmit 
7^0. 

The higher excitations in the two-electron quantum dots the lesser the 
influence of the Coulomb forces on the "crossing" of levels. For example, 
the value of the parameter k = l^/a* for singlet-triplet phase transitions 
decreases with increasing radial quantum number n. In particular, for the 
two-dimensional system we have obtained the following relation between pa- 
rameters k = lo/a* for a singlet-triplet transition at different n 

{h/a*)n=i ^ 2 + \m\ 
(/oK)„=o 7+|m| 

While the interplay between the magnetic field and the Coulomb forces de- 
termines the features of a phase transition (singlet triplet) for the ground 
state (n = 0) 0, mainly the magnetic field leads to the phase transitions for 
the high-lying states n > 0. 

The model allows the calculation of the magnetization M = —dE/dB. 
Since at low energy the magnetization is closely related to the slope of the 
ground energy, at T = K we obtain for n = 



dB 2 



x'^(\m\ + l) t 
m H — ^ ' 



xk t V{\m\ + \) J. / , g 



4 

(54) 



2v^(l + it2)3/4r(|^| + l)^'^^^^ + 2me^ ^ ' 

As it was mentioned in |^, |], the phase ground state transitions would be 
reflected in sharp discontinuities in the magnetization. The above exact 
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expression can be useful for the analysis of the experimental features related 
to the phase transitions. Also, it allows one to control the approximation 
made in the calculations in |TB . 



5 Summary 

Within the proposed model the analytical expressions for the energy levels 
and the magnetization of the two-electron quantum dots are obtained. The 
Coulomb interaction is treated exactly and from the analysis of the energy 
spectrum it follows that the interplay between the Coulomb forces and the 
magnetic field are an important ingredients for the prediction of the ground 
phase transitions. The Coulomb forces lead to the modification of the exter- 
nal potential and give rize to the effective confining potential of the relative 
motion. Their contribution in the properties of single-electron states de- 
creases with the increasing of the radial quantum number n. Finally, we 
would like to mention that the third extension (z) modifies the value of a 
magnetic filed needed to observe the phase transition: the thicker slab the 
larger value of a magnetic filed. We hope that the results presented here 
could be useful for the analysis of the electron properties in two-electron 
quantum dots and allow to make a conclusion on a deviation of the real 
confining potential from the harmonic oscillator one. 

A Two-dimensional case: matrix {n\hi\n) 

Here we describe some details of the calculations of the quantity 



oo 




< n\hi\n > 



I 



< n\ : 62 



2iT\/fujj(pqrj) 



n > (A.l) 
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Taking into account the following equations 



ika^+ —ika ;?+ 



e a ' e 



a + ik 



(A.2) 



the fact that 



[a a 



dr] 



-a(a+a+) 



a=0 



dot^ J \ x/K j 



(A.3) 



a=0 



and Eg. (pSf) , after some transformation we obtain 

, I -iB[a+^)-^Bia^) , , _ ^2 (^'^^)' + 4a/3)^' 

■ ■ '""^ "^"^ da-op- ■ j! ■ (1 - AaPY^'l' 

where B = t\/2. Using these results, we have for Eq. (|A.l| ) 



where 



and 



3 T{d/2 - 1/2) 



" 22p-'=r(A; + n-p + rf/2) 

Nk{n,d) = 



^o{n-p)\{2p-k)\{{k-p)\) 
In a particular case, n = 1, and n = 2 for S'^ we have 

4 



c _2 



^2 



ci(ci + 2) 



H 



a,/3=0 



(A.4) 
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B Three-dimensional case: definition of hj, 

Q(7) and Snh) 



Using the same technique as for the two-dimensional case and omitting te- 
dious calculations, we present the final result 



hi = hi + h2 + hs + hi 



■ ^2 



ho 



^ 1 + — ^ ■ ■ 

oo 



oo ■ — oo 



(1 + T^y/^ 



1 + 



r'^hu 



(1 + r 



2 

2\ ■ Pq ■ 



0r (1 + 7r2)3/2 (1 + r^)d/2+i " 



2 2 
Pq ■■ Qz 



hd — 



^^^^ -„2(l+r2)-r2(l+7T2) 



-2iTy/hw^{q^t) -2iTVnw(pqr)) 



(B.l) 



(B.2) 



(B.3) 
(B.4) 

(B.5) 



Q(7) = / 



0F(l + r2)'^/2 



.7r(|m| + l) 
2r(|m| + 1) ^ 



^n(7) = 



-2iT\/hjjj{pqri) 



< n\ : Cr, 



n > 



hT2) 



1 



(B.6) 
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7 r(n + i)r(| 


m\ 


+ 1) 


20F r(i + 


\m\ 


+ n) 



2n 

E( 

1=2 



-1)' 



r(i + o 



r(i + / + |m|) 



(B.7) 
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Figure Captions 

Fig.l The energy spectrum of a two-dimensional quantum dot in units 
of Hcuq as a function of the magnetic field strength uJc/oJq- The family of 
states with the quantum numbers A^ = 0, M = 0, n = and m < is 
shown (a) without and (b) including the Coulomb interaction between the 
two electrons. The arrow indicates the value of the magnetic field strength 
uJc/ojq = 1.91 where the second "crossing" occurs between the lowest states 
m = — 1 and m = —2. 

Fig. 2 Similar to the Fig.l for three-dimensional quantum dot (I/7 = 
Uz/ui = 3) including the Coulomb interaction between the two electrons. 
Here, the second "crossing" occurs at lJc/uj = 3.64. The decreasing of the 
ratio 1/7 leads to the " crossing" of levels at higher magnetic field strength . 
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